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Amplitude death is a dynamical phenomenon in which a network of oscillators settles to a stable 
state as a result of coupling. Here, we study amplitude death in a generalized model of delay- 
coupled delay oscillators. We derive analytical results for degree homogeneous networks that show 
that amplitude death is governed by certain eigenvalues of the network's adjacency matrix. In 
particular these results demonstrate that in delay-coupled delay oscillators amplitude death can 
occur for arbitrarily large coupling strength k. In this limit we find a region of amplitude death, 
which occurs already at small coupling delays that scale with 1/k. We show numerically that these 
results remain valid in random networks with heterogeneous degree distribution. 



INTRODUCTION 

Coupling individual oscillators can lead to amplitude 
death, the cessation of oscillations due to stabilization 
of a stationary state [1]. This phenomenon has been ob- 
served in chemical oscillators [2], electronic circuits [3], 
thermo-optical oscillators [4] , and coupled lasers [5] . Am- 
plitude death can be beneficial and is exploited for in- 
stance in feedback control applications [6], but can also 
be detrimental and even lethal if it occurred in interact- 
ing cardiac cells [7] . 

In mathematical models, amplitude death is commonly 
studied by investigating the conditions under which an 
unstable steady state of an isolated system is stabilized 
by coupling. It was found that for amplitude death to oc- 
cur, either the natural frequencies of the coupled systems 
need to be sufficiently disparate [8-11] or the coupling 
needs to be time-delayed [12]. 

While amplitude death can be observed already in 
systems of two oscillators, richer behavior is observed 
when complex networks of oscillators are considered. In 
such networks the nodes represent individual oscillators, 
whereas the links represent coupling terms. A central 
question is then how the dynamics are affected by net- 
work topology, the specific structure of nodes and links. 
For a given network of N oscillators this topology can 
be captured by the adjacency matrix A, a N x N ma- 
trix with Aij — 1 if a link exists from node j to node i 
and Aij — otherwise. A related matrix known to af- 
fect many network properties is the network Laplacian 
L = A — D. Here, D is a diagonal matrix, with Da = fej, 
the degree of node i, i.e. the number of nodes to which 
the focal node is connected. 

Delay induced amplitude death is often studied in sys- 
tems of coupled limit-cycle oscillators, such as Stuart- 
Landau oscillators. Studies on globally connected net- 
works [12] and rings [13] showed that amplitude death 
occurs inside islands in the parameter space of coupling 



strength and delay. While distributed delays seem to 
enlarge these islands [14], gradient instead of diffusive 
coupling has been found to suppress amplitude death in 
rings [15]. Further, transients between regimes of par- 
tial and complete amplitude death have been studied by 
numerical simulations [16]. Also the effect of delayed 
self-feedback has been studied in a single and two delay- 
coupled Stuart-Landau oscillators [17, 18]. 

Delay-induced amplitude death has been explored in 
detail in discrete-time maps [19-22] and in systems of os- 
cillators modeled by multi-dimensional ordinary differen- 
tial equations, such as the Rdssler and Lorenz oscillators 
[7, 23-25]. For a model of coupled maps, it was shown 
that the dynamics is governed by the largest eigenvalue of 
the network's Laplacian [19]. A similar result was shown 
in time-continuous system [7] . 

Previous studies have shown impressively that ana- 
lyzing amplitude death in simple systems, such as sim- 
ple oscillators or maps can yield a deep understanding, 
whereas coupling more complex systems such as Rossler 
or Lorenz oscillators points to additional complexities. A 
middle way is perhaps offered by using delay oscillators, 
such as Mackey-Glass [26, 27] and the Ikeda [28] oscil- 
lators. Because even a single- variable delay-differential 
equation (DDE) constitutes an infinite dimensional dy- 
namical system, a single DDEs can show sustained os- 
cillations, quasi-periodicity, and chaos. Studying ampli- 
tude death in delay-coupled systems of such delay oscil- 
lators offers the opportunity to consider a network of rel- 
atively complex coupled system, while keeping the equa- 
tions concise. 

Although a whole zoo of different synchronization 
types have been found for systems of coupled delay- 
oscillators without coupling delay [29-34], delay-coupled 
delay oscillators, have been less explored. One excep- 
tion is [6] where a system of two mutually delay-coupled 
delay oscillators was investigated numerically, but more 
complex coupling topologies were not considered. 
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A model for complex networks of general scalar delay- 
coupled delay oscillators was proposed in two previous 
publications of the present authors [35]. The latter pa- 
per highlighted that certain meso-scale structures in net- 
works have a distinct impact on the network-level dynam- 
ics, whereas the former discussed the effect of the degree 
distribution, the probability distribution of node degrees. 
For degree- homogeneous networks (DHONs), where ev- 
ery node is connected to the same number of of other 
nodes, the dynamics of delay-coupled delay oscillators 
were studied analytically. This revealed the bifurcation 
lines at which the dynamics of the system change qual- 
itatively, and linked them to certain eigenvalues of the 
networks adjacency matrix. For degree-heterogeneous 
networks (DHENs) , numerical explorations showed qual- 
itatively similar dynamics, with minor corrections. 

In the present paper, we use the previously proposed 
model to study amplitude death. For this purpose, we 
focus on a regime where the uncoupled systems show 
non-stationary dynamics. Then, we study the parameter 
space of coupling strength and coupling delays to find 
stable areas, which correspond to regions of amplitude 
death. For DHONs, we find that there exists a region of 
amplitude death for coupling delays larger than a certain 
threshold value, which approaches zero as l/(dk) where d 
is the common degree of nodes and the coupling strength 
k approaches infinity. By using a numerical sampling 
method, we confirm this result for degree-heterogeneous 
networks, where d is now the mean degree of the network. 
This shows that in the limit of large coupling strength the 
dynamics is governed by a global property of the network, 
while all other topological properties become negligible. 

MODEL 

We consider networks of N nodes i, carrying a dynam- 
ical state X iy which can represent for instance an ecolog- 
ical population or the abundance of RNA molecules in 
a genetic oscillator. The variable is subject to gain and 
loss terms, where the loss is instantaneous, while the gain 
term of Xi at time t depends on the value of Xi at time 
t — t. Furthermore, X^ continuously diffuses to topolog- 
ical neighbors causing an additional loss. The exported 
material arrives at the respective neighboring node j af- 
ter a travel time delay (5. Thus, the dynamics is governed 
by 

Xi = G(XT) - L( Xl ) + {A l3 F{X s 3 ) - AjiFiXi)) , 

3 

(1) 

where G, L and F are positive functions describing 
growth, loss, and coupling, respectively. Instead of re- 
stricting these functions to specific functional forms, we 
consider a general model comprising the whole class of 
models that includes well-studied examples such as the 



delay-coupled Mackey- Glass and the Ikeda systems. 

Here we only consider networks for which the number 
of outgoing links d° ut equals the number of incoming links 
for each node i, such that d° ut = d° ut =: di, where di 
is the degree of node i. We note that this class includes 
among all bidircctionally coupled networks. For all net- 
works in the class considered, a homogeneous steady state 
Xi* = X* exists, for each steady state X* of the isolated 
system, which satisfies G(X*) = L(X*). 

We apply the method of generalized modeling [36, 37], 
which analyzes the dynamics of general system by a 
direct parametrization of the Jacobian matrix J, with 
Jij = dXi/dXj in X* . This matrix constitutes a local 
linearisation and thus governs the dynamics of the sys- 
tem close to a steady state under consideration [38]. In 
order to express the Jacobian in terms of interpretable 
parameters, we normalize the system to an arbitrary 
homogeneous steady state Xi* and introduce normal- 
ized variables x\ = Xi/X* and normalized functions 
f(xi) — F(xiX*)/ F(X*), denoted by lower case sym- 
bols. Using X* = X T * = X s *, we rewrite Eq. (1) as 

±i = a(g(xT) - l( Xi )) - dip/fa) +^^^ J /(x«). (2) 

3 

The quantities a = G(X*)/X* = L(X*)/X* and /3 = 
F(X*)/X* are unknown constants and can thus be in- 
terpreted as parameters of the Jacobian. These parame- 
ters have the unit of inverse time and describe turnover 
rates. In the following we set a = 1 by a time scale 
normalization. 

In order to analyze the stability of the homogeneous 
steady state, we proceed as in [39] and assume that, close 
to the steady state, the DDE has exponential solutions 
y(t) = y(0) exp(At), where y(t) = x(t) — 1 is a small per- 
turbation from the steady state at x = 1. Using this 
ansatz together with the linearization of Eq. (2) the Ja- 
cobian matrix has the entries 

Jf := g' exp(-Ar) - /' - di/3f for i = j 
J° := ' Pfexp(-\5) tovAij = l (3) 

otherwise 

where A is an eigenvalue of the Jacobian and thus a solu- 
tion of the characteristic polynomial P(X) = det(AI— J). 
The quantities g' , I', f denote derivatives of the normal- 
ized functions in the steady state, or equivalently log- 
arithmic derivatives of the original functions, i.e. g' = 
dg/dx\ x = Slog (G)/d\og(X)\^. Although we do not 
specify the functions or the steady state under consider- 
ation, these derivatives are formally constants and can 
thus be interpreted as unknown parameters of the sys- 
tem. Closer inspection reveals that these so-called elas- 
ticities are generally easily interpretable in the context of 
applications and have a number of additional benefits [36] . 
Since /' only appears in products together with /3, we in- 
troduce the effective coupling strength k :— /3f. 
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The steady state under consideration is asymptotically 
stable if all roots of the characteristic polynomial have 
negative real-parts and it is unstable if at least one root 
has a positive real part [1]. In DDEs the explicit ap- 
pearance of A in the Jacobian turns the characteristic 
polynomial into a transcendental function, which may 
have an infinite number of solutions. To make progress 
analytically it is therefore advantageous to directly com- 
pute the bifurcation points in parameters space at which 
eigenvalues acquire positive real parts. This is general 
simpler than computing the spectrum for a single pa- 
rameter set. In numerical explorations the computation 
of the spectrum can be avoided by checking for positive 
eigenvalues with a method based on Cauchy's argument 
principle [35, 40]. 



BIFURCATIONS OF DEGREE HOMOGENEOUS 
NETWORKS 

In a previous publication [35], we demonstrated that 
the characteristic equation P(X) — of DHONs can be 
decomposed into N independent equations, each of which 
corresponds to a single eigenvalue of the adjacency ma- 
trix. In this section we briefly recapitulate this derivation 
that will form the basis for our discussion of amplitude 
death below. Because the treatment involves both the 
eigenvalues Aj of the system's Jacobian matrix, and the 
eigenvalues Ci of the system's adjacency matrix, we re- 
fer to Aj as dynamic eigenvalues and to Cj as topologial 
eigenvalues. 

In a DHON all diagonal elements Jf are identical. This 
allows us to substitute A — J d in the characteristic poly- 
nomial P(A) by z, which yields P'{z) = det(zl — J°A). 
The roots of P'(z) are given by Cj J°, where Cj denotes one 
of the N topological eigenvalues. The back-substitution 
yields N independent equations, 



A = J d (A) +qJ°(A). 



(4) 



We note that the topological eigenvalues c are in general 
complex, so that c = |c|e 1 ^ , where ip c is the complex 
phase of the topological eigenvalue. 

For determining the stability of the system, we recall 
that the stability changes at bifurcation points where a 
dynamical eigenvalue crosses the imaginary axis and be- 
comes purely imaginary, so that A — iuj. By separating 
the real and imaginary part of Eq. (4) we obtain 

= g' cosO) -l'-dk+ |c|fccos(», (5) 
u = — </sin(</>) — |c|fcsin(^), (6) 

with <p = ujt and ip := u)5 — ?p c . 

The equation system contains three unknown vari- 
ables, <f>, ip and ui. Given a solution triplet (<p,ip,u>), we 
can find another solution {—<p,—ij),—u}). The solutions 
with negative oj trace the same bifurcation lines as those 



with positive us. We can thus restrict the analysis to so- 
lutions uj > without missing bifurcations. Further, we 
are only interested in solutions with r > 0, and therefore 
(f) > 0. By choosing <f> as a free parameter, we can find 
a parametric representations of the bifurcation lines, but 
need to distinguish three cases depending on the value of 
the topological eigenvalue c. 

For the case c = the Eqs. (5,6) are independent of 
ip and 5, respectively. Thus, the bifurcations are vertical 
lines in the (fc, 5)-plane, with 



1 h 



(7) 



where h := (</ cos(0) — V)t and (f> has to satisfy f((j>) := 
4> + g'r sin((/>) = 0. The equation f(4>) = may generally 
have several solutions. However, for the parameters used 
throughout this paper the solution is unique. 
For the case |c| = d, wc find 



T 



2dh 



(8) 



In order to obtain positive solutions, h((f>) needs to be 
positive. Therefore, we have to restrict <f> to the intervals 
1* = with 0* = 2irr+(f> k , = 2vr(r+l)-0 fe and 

(\> k = cos -1 (/'/(/). The parameter r that appears here is 
the first of two integer parameters that we introduce to 
enumerate the bifurcation lines. 

Finally, for the case < |c| < d, wc find 



ki,2 = 



ldh±y / 



.2 ' 



with 



a = d 2 h 2 -(d 2 -\c\ 2 )(f + h 2 ). 



(9) 



(10) 



For obtaining real and positive solutions, we require not 
only h((f>) > but also a(<f>) > 0. This imposes additional 
constraints on (f> that can be calculated numerically. In 
contrast to the case |c| = d, valid solutions can only be 
found inside a finite number of intervals If, 

We remark that the case |c| > d does not need to 
be considered because the eigenvalues c of the adjacency 
matrix cannot exceed the highest node degree, which is 
d in DHONs. 

For the cases above, we compute the corresponding val- 
ues of S at which the bifurcation occurs from Eqs. (5,6). 
We find 



rL,R 



L,R 



± COS 



h 



(ii) 



(12) 



where s, is the second integer parameter introduced to 
enumerate solution branches. Furthermore, the indices 
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L, respectively R, are introduced to denote the positive, 
respectively negative, sign branch of ip. To obtain phys- 
ical solutions we consider the L-branch for all f(<p) < 
and the R-branch otherwise. Evaluating these, we find 
positive values of S for non-negative integers s. 

One implication of Eq. (12) is revealed when we con- 
sider that the arcus cosine needs to be smaller or equal 
to one. This is only possible if k < k max with 



pan _ _ 



1 + 9 



(13) 



lines corresponding to topological 
cj < d can only be found for finite 



Thus, bifurcation 
eigenvalues c with 
values of k. Hence, only eigenvalues with \c\ — d can af- 
fect the stability for sufficiently large coupling strength, 
independently of the coupling topology. 



AMPLITUDE DEATH IN DEGREE 
HOMOGENEOUS NETWORKS 

Having derived the results in the previous section, we 
now turn to the analysis of amplitude death. In the pa- 
rameter range g' < — V < 0, we find that an isolated node 
is unstable if t > r* where r* = r* n and 



2nr 



(14) 



The isolated system thus exhibits non-stationary (e.g. os- 
cillatory) dynamics when the reproductive delay r is cho- 
sen sufficiently large. The stationary solutions can then 
potentially be stabilized by coupling the oscillators. In 
general stabilization will depend on the coupling topol- 
ogy, coupling delays <5, and the coupling strength k. In 
the following we seek to identify the effect of the cou- 
pling topology on the areas in the (fc, 5)-space where the 
homogeneous steady state is stable, such that amplitude 
death can occur. 

Our analytical treatment above has identified a con- 
nection between the topological eigenvalues and the dy- 
namical stability. For illustration of the analytical results 
we consider two specific topologies: A fully connected 
networks of 3 nodes and a ring of four nodes. 

Let us start with the fully connected network (Fig. la). 
All such networks have a topological eigenvalue c = d = 
N — 1 and a d-fold degenerate topological eigenvalue c = 
— 1. For c = — 1 computation of the bifurcation lines only 
yields physical solutions for the r = branch. In this 
branch different values of s generate different segments 
of bifurcation lines, which connect such that a single long 
bifurcation line running from S = to S = oo is formed. 
Because |c| < d, all bifurcation points in this line occur 
at finite values of the effective coupling strength k. 

For the remaining eigenvalue c = d, the situation is 
more complex. Here, computation of the bifurcation 



M 

o 
o 



o 




o 



\ 2 

.2 

M 













/ 1 — 














1 1 ' ' — 

i 


1 1 





12 3 4 
coupling strength k 



-*;:ttitiitttttitiifl 



0.5 1 10 20 

elasticity ratio —l/g tongue index r 

FIG. 1. Bifurcation diagrams indicating the dynamics for a 
fully connected network of 3 nodes (a) and a ring of four nodes 
(b). Lines mark bifurcation points corresponding to different 
topological eigenvalues Ci (see text). Labels of the form [r, s) 
on some of the lines state the corresponing indices r and s 
enumerating the analytical solution branches. For the limit 
of high coupling strength the bifurcation lines asymptotically 
approach limiting values, some of which are indicated on the 
right border of a) by symbols of the form <$J?'* , where 6 and 
t denote bottom and top branches, respectively. Gray dots 
mark stable states that were found numerically by uniformly 
sampling the (k,S) plane. The two lower plots (c,d) show 
the dependence of limiting values <5j?,'s on the ratio of local 
parameters — l/g (c) and the dependence of <5j? r+1 , <5* r on the 
index r. Parameters: g' = —1,1' = 0,r = 5. 



points yields a family of separate tongues corresponding 
to different values of r and s. Similar families of bifurca- 
tion lines have been observed in other delay systems [41]. 
In the present paper, each of the tongues consists of two 
segments, to which we refer as the bottom and top part 
respectively. An exceptional solution is again the case 
r = 0. Here, the lower boundary of the tongue is formed 
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by the top branch of the (r = 0, s = 0) solution, while the 
upper boundary is formed by the bottom branch of the 
(0, 1) solution. One can thus think of the r = case as 
an inside-out tongue, which provides the right intuition 
for understanding the results below. 

Let us now turn to the ring network (Fig. 1(b)). All 
even rings are bipartite networks, which means that it is 
possible to color the nodes in two colors such that no node 
has a topological neighbor of the same color. Bipartite 
degree homogeneous networks have topological eigenval- 
ues at c = d and c = — d, additionally the 4-ring has a 
two-fold degenerate eigenvalue at c = 0, which results 
from symmetry. For the topological eigenvalue c = 
only the (r = 0,s = 0) solution corresponds to a physical 
bifurcation line. This line occurs at a constant value of 
the coupling strength. The eigenvalues with c = d and 
c = — d each generate a family of tongues that are sim- 
ilar to those in the fully connected network. However, 
the tongues with c = — d are shifted such that they are 
centered on the gaps between the tongues with c = d. 
From the corresponding eigenvectors one can see that 
the family of tongues with c — d corresponds to an insta- 
bility with respect to in-phase oscillations, whereas the 
tongues with c = —d correspond to anti-phase oscilla- 
tions. It is intuitive that this instability arises from the 
eigenvalue c = — d which is directly linked to bipartite- 
ness, because none-bipartite networks could not sustain 
anti-phase oscillations. 

The analytical results above can be confirmed by a 
numerical sampling procedure, in which we pick param- 
eter sets uniformly from the (fc,<5)-plane and evaluate 
their stability numerically [35]. For visualization we plot 
only those points that are found to be stable. The re- 
sults shown in Fig. 1 shows that such stable parameter 
combinations are found only in certain regions that are 
sharply delineated by the theoretically predicted bifurca- 
tion lines. The figure shows that stability requires that 
the respective point lies inside all tongues with r = 
but outside all tongues with r ^ 0, which confirms the 
inside-out nature of the tongues with r = mentioned 
above. Thus for stable points, the bifurcation lines corre- 
sponding to r = impose a minimum coupling strength, 
whereas the solutions with r =^ impose a maximum 
coupling strength. 

We note that the c = topological eigenvalue never 
causes a bifurcation of a stable solution. This eigen- 
value gives rise to a physical branch with r — and 
thus potentially imposes a lower limit for the coupling 
strength in stable states. However, this limit is always 
below a higher limit imposed by the r = branches of 
other eigenvalues, and thus has no direct relevance for 
stability. This conforms to the general observation that 
only the largest positive and the smallest negative topo- 
logical eigenvalue cause bifurcations that border stable 
regions. For DHONs it is probably possible to prove this 
rigorously, but the proof is beyond the scope of the cur- 



rent paper. Nevertheless, it is interesting to note that 
in all DHONs the largest positive eigenvalue is c = d, 
whereas the smallest negative eigenvalue is maximal for 
fully connected networks (c = 0) and minimal for bipar- 
tite networks (c = —d). This illustrates why amplitude 
death is most likely in the fully connected network and 
least likely in bipartite networks. All other networks fall 
between these two extreme cases, which motivated our 
choice of examples. 

STRONG COUPLING LIMIT 

In the following, we focus on large coupling strength 
k. Above we showed already that only those bifurca- 
tion lines that correspond to topological eigenvalues with 
|c| = d can extend to arbitrarily large coupling strength. 
Except in bipartite networks, this condition is only sat- 
isfied for the eigenvalue c = d. In the following we thus 
focus solely on this eigenvalue. 

In the limit k — ¥ oo, the bifurcation lines approach 
constant values of the coupling delay 5}?'*, respectively, 
where t and b denote values for top and bottom branches. 
For calculating 6~'g, we study Eq. (8) and find that k 

becomes infinitely large if <j) approaches and Using 
Eq. (11) we obtain 

t 2irs + r b 2ns + 

r ' s 2nr + <t> k l r ' s 2n(r + 1) - <p k 1 ' 

We note that both equations scale linearly with r, such 
that all tongues scale with the internal delay of the os- 
cillators. Furthermore, recalling that <jr = cos _1 (Z' /</), 
we see that for g' = —I', the bottom and top limits are 
identical for each bifurcation line {r,s), such that the 
unstable areas disappear (see also Fig. 1(c) for selected 
lines). 

Let us now investigate the question whether in the 
limit of large coupling strength amplitude death is still 
possible. Already visual inspection of Fig. 1(a) reveals 
that there are large channels of stability between the 
tongues, which seem to extend to high values of k. In 
particular, notable is the l:l-resonant channel around 
5 = t and the bottom channel at small values of S. As 
we increase the coupling strength, all stable channels are 
successively narrowed down as new tongues, correspond- 
ing to branches with higher r, become relevant. For in- 
stance, the l:l-resonant channel is bordered from below 
by the top branch of (r, r) and from above by the bottom 
branch of (r, r + 1); the bottom channel is bordered from 
below by the top branch of (0, 0) and from above by the 
bottom branch of (r, 1). 

We now focus specifically on the bottom channel, 
which comes arbitrarily close to the case of zero cou- 
pling delay and thus the well studied case of networks 
of delay oscillators without coupling delay. We start by 
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FIG. 2. Amplitude death for large coupling strength, (a) 
Bifurcation lines from the eigenvalue c — d (solid lines) and 
c = —d (dashed lines). The approximate tip positions of 
tongues have been marked by circles (c = d) and squares 
(c = — d). The orange dashed line marks the analytical ap- 
proximation to the bottom boundary of the stable region. 
For sufficiently large k, all degree homogeneous networks are 
stable between this bifurcation line and the tongues corre- 
sponding to c = — d. Non-bipartite networks are even stable 
up the tongues corresponding to c = d. (b) Numerical results 
from uniform sampling of the space shown for random net- 
works with TV = 10 nodes and K = 15 edges. Red dots mark 
unstable parameters points, such that the large white regions 
mark stable areas that correspond to amplitude death. In 
particular, for sufficiently k, no unstable networks are found 
for 5 between the bottom bifurcation line with (d = 2K/N) 
and the tip positions for c = — d. Most networks are stable 
up to the tip positions for c = d. In the double logarithmic 
plot the region of amplitude death growth with increasing k 
as its boundaries approach power-laws with exponent —1 and 
— 1/2 respectively, in agreement with analytical predictions. 
Parameters are g' = — 1, 1' = 0,t = 5. 



considering the lower boundary of this channel by study- 
ing the asymptotic behavior of the (0, 0) bifurcation line. 
We note that in the relevant limit this is given by the 



L-branch of Eqs. (11,12) with 4> approaching k . In this 
case h — > 0, so that we can approximate Eq. (12) by us- 
ing cos _1 (l 



x) « -v/%1- With T * = <t> 7vV 2 - l' 2 and 
k = cos _1 (/'/(7'), we find the critical coupling delay 



l 

dk' 



(16) 



where we dropped the indices for simplicity. This result 
shows that the critical S, which marks the lower boundary 
of the bottom channel, is inversely proportional to the 
effective coupling strength dk. 

For finding the critical 5 that marks the upper bound- 
ary of the bottom channel, we approximate the tip posi- 
tions of the tongues (r, 1) in the limit of large r. In this 
case we can approximate Eq. (8) with 



2 
2dhr' 



(17) 



Using that <p is large, we find local minima at (f> = 
(2r + l)7r. Inserting this result into the R-branch of 
Eqs. (11,12) yields an expression for the 8- value. For 
large k, we see from Eq. (11) that 6 is inversely propor- 
tional to <f>, while we see from Eq. (17) that k is propor- 
tional to (f> 2 . Thus, 6 oc 1/Vk. 

In summary, these results show that in the limit of 
large coupling strength there is a region of amplitude 
deaths at very small coupling delays. The onset of am- 
plitude death occurs at a value of the coupling delay that 
scales as 1/k, whereas the width of the amplitude death 
region (in terms of the coupling delay) scales as \j\fk. 
We note that the reasoning presented in this section is 
not limited to the example topologies discussed above, 
but holds for all degree-homogeneous networks. While 
additional topological eigenvalues with |c| = d could ex- 
ist in directed or bipartite networks, the corresponding 
tongues would narrow the channels of amplitude death 
by a factor, but would show the same scaling behavior. 



AMPLITUDE DEATH IN HETEROGENEOUS 
NETWORK 

As the final step of our analysis we investigate ampli- 
tude death in degree heterogeneous networks (DHENs). 
In particular we focus again on the limit of large coupling 
strength and small coupling delays that corresponds to 
the bottom channel computed analytically in DHONs. 
For DHENs comparable analytical calculations are not 
easily possible. We therefore explore the stability of these 
systems numerically. Figure Fig. 2 shows a scatter plot in 
which each point correpsonds to specific values of k and 
S, sampled uniformly, and a random Erdos-Renyf net- 
work with N = 10 nodes and K = 15 links. For each of 
these randomly generated sample networks we used the 
numerical method described in [35, 40] to compute the 
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coupling strength k 



FIG. 3. Onset of amplitude death. Scatter plot of parameter 
values for which one randomly drawn network with N = 10 
nodes and K links is stable and another one is unstable. 
Thus the dots mark the border between stable and unstable 
regions. The borders of network ensembles with K = 14, 
K — 15 and K = 16 approach distinct lines, which are 
given by Eq. (16) with d — 2K/N (horizontal lines). The 
S- values are normalized by the line for d = 3. Parameters are 
g' = -l,l' =0,t = 5. 

stability. In the figure unstable samples where marked 
by a red dot, whereas stable regions remain white. 

From the scatter plot we see that for sufficiently large 
k, no unstable networks can be found in a large region at 
high coupling strength and small coupling delays. This 
region corresponds remarkably well to the analytical re- 
sults that have been obtained for DHONs. In particular 
there is a sharp loss of stability at he bottom bifurca- 
tion line. For sufficiently high coupling strength almost 
all unstable samples fall into the tongues of instability 
generated by the eigenvalue c = d in DHONs, whereas 
few fall into the tongues with c = — d. These samples 
can probably be explained by the random creation of bi- 
partite networks, and would hence be virtually absent in 
larger networks. 

As a further test we study the onset of amplitude death 
as the bottom (0,0) bifurcation line in more detail. To 
capture the onset of amplitude death we are interested 
in those systems values where for given values of 5 and 
fc the stationary solution is stable in one network topol- 
ogy and unstable in another topology. We thus repeat 
the sampling procedure above, but draw two topologies 
for every 5, k value pair. We discard all samples except 
those where the stationary state is stable in one topology 
and unstable in the other. Since we already know that 
the stability transition at the bottom bifurcation line is 
narrow for large k we sample this region selectively by 
drawing 8 from a bounded uniform distribution centered 
on (5(2=3 (fc), the bifurcation point for DHONs with degree 



d — 3, which is given by Eq. (16). This results in a uni- 
form sampling of the parameter plane shown in Fig. 3. 
We observe that for large k the transitions occur only at 
specific values of 6. These values very closely approxi- 
mate the analytical predictions of the bifurcation points 
in DHONs, where we replaced the homogeneous degree 
by the mean degree d = 2K/N of the respective DHEN. 
For large coupling strength the bifurcation line that gov- 
erns the onset of amplitude death at small coupling de- 
lays thus seems to depend only on the mean degree of 
the network, whereas all other topological properties, at 
least in the Erdos-Renyi ensemble, can be neglected. 

CONCLUSIONS 

In the present paper we investigated amplitude death 
in general networks for delay-coupled delay oscillators. 
Building on a previous result we were able to study the 
regions for amplitude death analytically for degree homo- 
geneous networks. In particular we considered the limit 
of large coupling strength for which we showed that all 
relevant bifurcation lines can be traced back to certain 
eigenvalues of the systems adjacency matrix. 

We showed analytically that in degree homogeneous 
networks regions of amplitude death exist at large cou- 
pling strength fc. Specifically, we investigated a region 
for which the onset of amplitude death occurs already at 
very small values of the coupling delay that scales as 1/fc. 
The width of this region scales as 1 j\fk and thus is signifi- 
cant, at least in a logarithmic sense. It is remarkable that 
this region of amplitude death comes arbitrarily close to, 
but never reaches the case of undelayed coupling. 

The results obtained here for delay oscillators contrast 
with the well-studied case of Stuart-Landau oscillators 
[12], which do not show amplitude death for arbitrary 
large coupling strength and arbitrary small coupling de- 
lays. We can speculate that the strong self-feedback in 
delay oscillators might be in destructive resonance with 
the terms arising from the coupling to neighboring nodes. 
Thus, in contrast to Stuart-Landau oscillators a direct 
force working against the oscillatory dynamics might not 
be necessary. 

Numerical investigations indicate that the analytical 
results for degree homogeneous networks can be extended 
to the case of degree heterogeneous networks. For suffi- 
ciently large coupling strength we found that numerical 
results for degree heterogeneous networks were in per- 
fect agreement with analytical expectations for degree 
homogeneous networks. While these results were based 
on network topologies from an ensemble of Erdos-Renyi 
random graphs one can assume that they should hold 
generally, at least for networks with exponentially decay- 
ing degree distribution. 

Regarding degree heterogeneous networks we noted 
that the onset of amplitude death for strong coupling 



and small coupling delays seems to depend on the net- 
works mean degree, but not on other properties. This 
is remarkable as it shows that the node dynamics be- 
come sensitive to a global and hence delocalized quantity, 
which hints at a diverging correlation length. This can be 
made plausible by considering that strong coupling cre- 
ates a stiff system that can rapidly communicate over the 
relatively small diameter of the networks. Furthermore, 
we have seen in the degree homogeneous networks that 
the dominating instability in the limit of strong coupling 
are in-phase oscillations, which are an inherently global 
phenomenon. 
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